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Abstract 

The dynamics of entanglement and fidelity for a subsystem of two separate spin-1/2 qubits 
prepared in Bell states is investigated. One of the subsystem qubit labelled A is under the influence 
of a Heisenberg XY spin-bath, while another one labelled B is uncoupled with that. We discuss 
two cases: (i) the number of bath spins N — > oo; (ii) N is finite: N = 40. In both cases, the 
bath is initially prepared in a thermal equilibrium state. It is shown that the time dependence 
of the concurrence and the fidelity of the two subsystem qubits can be controlled by tuning the 
parameters of the spin-bath, such as the anisotropic parameter, the temperature and the coupling 
strength with qubit A. It is interesting to find the dynamics of the concurrence is independent of 
four different initial Bell states and that of the fidelity is divided into two groups. 
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I. INTRODUCTION 



Entanglement, which exhibits a very pecuhar correlation among the degrees of freedom 
of a single particle or the distinct parts of a composite system, is the most intriguing 
feature of quantum composite system and a vital resource for quantum computation and 



quantum communication 



ite system ai 

2l,SQ,y,B- 



In the field of quantum information theory, it is a 



fundamental issue to create, quantify, control and manipulate the entangled quantum bits 



which are often composed of spin-half atoms in different problems 



lots of works 



Id, 



11 



Particularly, 



12l | are devoted to steering two initially entangled qubits through 
an auxiliary particle or field (for instance, another spin qubit or a bosonic mode), which 
interacts with only one of them. It is a very exciting motivation, yet those approximated 
models neglect the actual complex environment of the quantum qubits. And it remains 
an important open question how the entanglement degree responds to the influence of 



environmental noise 13|]. Practically, the spin qubits are indeed open quantum subsystems 
and exposed to the influence of their environments Q, H, 0, Il7l |. In most conditions, the 
coupling between the subsystem and environment will degrade the entanglement degree 
between the subsystem qubits. In some other conditions, however, a specially structured 
and well designed bath can be conceived as a protection device to suppress the negative 



influence from itself or other noise sources 



18 



19 



20| 



For the spin subsystem, there are two important modes of bath: (i) boson-bath, e.g.. 



the Caldeira-Leggett model 



21| : (ii) spin-bath, e.g., the model used in Ref. [22|. Here we 



discuss the latter one. It is well known that the localized spins in solid state nano-devices, 
;he most promising candidates for qubits due to their easy scalability and controllability 
3, S]> mainly subject to the influence from the nuclear spins, which constitute a type 
of spins-1/2 environment. It is an almost intractable computation task to deal with such 
a spin-spin-bath model for its giant number of degrees of freedom. Therefore, physicists 
resort to some approximations or simplifications, such as the Markovian 2^ schemes and 



the non-Markovian ones 



25(1, which have been developed in the past two decades. Based 



on these schemes, plenty of analytical and numerical methods were exploited to study the 
reduced dynamics of subsystem consisted of spins-1/2 by tracing out the degrees of freedom 
of the spin-bath. Some recent works focused on the center spins in a network configuration. 
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in which the form of bath is specially structured, such as a thermal bath 26|, |27|] , a bath 
via Heisenberg XX couplings [28| and a thermal spin bath via Heisenberg XY couplings 

3,Q- 



In this paper, an open two-spin-qubit subsystem (two qubits labelled A and B respec- 
tively) is explored as a target quantum information device with a spin bath of a star-like 



17 



But there are 



configuration. The model is something like the one considered in Ref. 
significant differences between them. It is supposed that at the beginning, the subsystem is 
prepared as one of the Bell states (Einstein- Rosen- Podolsky pairs) [3ll, l32] : 



|ei) = l/v^(|ll) + |00)) 

les) = l/v^(|10) + |01)) 

lea) = l/v^(|ll)-|00)) 

|e4) = l/v^(|10)-|01)) 

Then Qubit B is moved far away or isolated so that not only the coupling between Qubit 
B and A, but also the interaction of B with the spin-bath could be neglected. The bath is 
regarded as an adjustable auxiliary device to control the time evolution of the two qubits. 
And the evolution is represented by the concurrence and fidelity dynamics as functions of 
various parameters associated with the bath. The reduced dynamics of the subsystem is 
obtained by a numerical scheme combined by the Holstein-Primakoff transformation and 
the Laguerre polynomial expansion algorithm. We consider two conditions, in which the 
number of spins in the bath is infinite and finite. The rest of this paper is organized as 
following. In Sec. [Ill we first give the model Hamiltonian and its analytical derivation; and 
then we introduce the numerical calculation about the evolution of the reduced matrix for 
the subsystem. Detailed results and discussions are in Sec. IIIII And conclusion is given in 
Sec. Ml 
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II. MODEL AND METHOD 



A. Hamiltonian 

The subsystem is consisted of two entangled spin-1/2 atoms labelled A and B respectively, 
between which there is no coupling. Qubit A interacts with a spin-1/2 bath via a Heisenberg 
XY interaction while B does not. The total system Hamiltonian, similar to those considered 



in Refs. 



17 



28|, 



33|, is divided into Hs, Hb and Hgb- They represent the subsystem, the 



bath and the interaction 



3, 



17 



3J] part between the former two terms respectively: 



H. 



sb 



H, 



9o 



N 



2VN 



i=l 



9 

2N 



N 



(1) 
(2) 

(3) 



where /xq is half of the energy bias for the two-level atom A or B; go is the coupling strength 
between the subsystem spin A and the bath spins; g represents the mutual interactions 
among the bath spins. 7 (0 < 7 < l)_is the anisotropic parameter. When 7 = 0, the XY 



interaction is reduced to an XX one 
the Pauli matrices. 



17|. The X, y and z components of the matrix a are 




a" 



-i 

1 



1 
-1 



(4) 



The indices i in the summation run from 1 to iV and is the number of the bath spins. 



Adopting (7^ = l/2(cr^' ± ia^), we can rewrite the Hamiltonians ([2]) and ([3]) by: 



H. 



sb 



r N 



N 



i=l 



N 



(5) 



(6) 



Substituting the collective angular momentum operators J± = J2iLi '^t |25| into Eqs. ([5 



and (Q, we get 



Hsb 



9o 



^ [7 ( J+ J+ + J-J-) + {.hJ- + J-J+- N)] . 



By the Holstein-Primakoff transformation 361], 

J+ = b+{VN-b+b), J_ = (ViV - 6+6)6, 

with [6, 6+] = 1 and the first order approximation of 

'N-b+b . b+b 



(7) 
(8) 

(9) 



1 - 



1 - 



n 



the Hamiltonians Eq. 

Hsb = go 



N 2N 2N' 

and Eq. (IHl) can be written as: 



(b' 



--9\ 7 



{b' 



)(6+ ) + {b )ib 



(10) 



+ (b^ )(b ) + (b )(b^ 



--9\ 7 



+ 2h + 



n{2n^ 



%Nh - n + 1) 



b'^h 

n-2 
2N 



h — 1, 
2N ' 



(11) 



4Ar2 



Utihzing the collective environment pseudospin J and the Holstein-Primakoff transfor- 
mation, one could reduce a high-symmetric spin bath, such as the one we considered, into 



a single-mode bosonic bath field ly, [17(1 . The transformed Hamiltonian is just like a spin- 



boson model in the field of cavity quantum electrodynamics (CQED). And the effect of the 
single-mode bath on the dynamics of the two subsystem qubits is interesting although the 
bath only directly interacts with one of them. The model might be helpful to understand 
the magic essence of quantum entanglement and practical in manipulating the quantum 
communication. 



B. Calculation method 



The whole state of the total system is assumed to be separable before t = 0, i.e. p(0) = 
[■^(O)) (?/'(0)| ® Pb- The subsystem [^/^(O)) is prepared as one of the four Bell states, \ei),i = 

5 



1,2,3,4. The bath is in a thermal equihbrium state, Pb{0) = e~^''^^^'^ /Z , where Z = 
Tr (^e"^''/'^'^^) is the partition function. The Bohzmann constant fc^ is set to be 1 for the 
sake of simphcity in later calculation. To derive the density matrix p{t) of the whole system, 

p{t) = exp{-iHt)p{0) exp{iHt), (12) 

we need to consider two factors. 



i) To express the thermal bath, we use the method suggested by Tessieri and Wilkie 



18|,l20|,l29|: 



N 



Ph 



(0) = ^ |0m)c^m(0m|, 



m=l 

p — Em /T 



N 
m=l 



(13) 

(14) 
(15) 



where \4>rn), m = 1,2,3, N, are the eigenstates of the environment Hamiltonian Hh, and 
Em the corresponding eigen energies in increasing order. On the condition of thermodynam- 
ics limit, i.e. oo, Eq. (fTOl) and Eq. (ITT!) are simplified as: 



Hb = g[iib+^ + b^) + 2b+b]. 



(16) 
(17) 



Then in Eq. ( |T3i) and Eq. ( |T5i) should be replaced with a cutoff M linking to a certain 
high energy level. By the above expansion, the initial state can be represented by: 



p(o) = 5^^^i*^(o))(*^(o)i, 1*^(0)) = \^m\K)- 



(18) 



(ii) For the evaluation of the evolution operator U (t) 



e xp 



iHt) , we apply the Laguerre 



polynomial expansion scheme, which is proposed by us 20|, |29|, |37| , into the computation. 

a+l oo / .j^ ^ A; 



U{t) 



1 + it 



k=0 



it 



+ it 



(19) 



L^{H) is one type of Laguerre polynomials |38| as a function of H, where a (—1 < a < oo) 
distinguishes different types of the Laguerre polynomials and k is the order of them. In 
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real calculations the expansion has to be cut at some value of fcmax, which was optimized 
to be 20 in this study (We have to test out a fcmax for the compromise of the numerical 
stability in the recurrence of the Laguerre polynomial and the speed of calculation). With 
the largest order of the expansion fixed, the time step t is restricted to some value in 
order to get accurate results of the evolution operator. At every time step, the accuracy 
of the results will be confirmed by the test of the numerical stability — whether the trace 
of the density matrix is 1 with error less than 10~^^. For longer time, the evolution can 
be achieved by more steps. The action of the Laguerre polynomial of Hamiltonian to the 

aguerre polynomial. The scheme is of an 



39 



40[ | and is pretty well suited to many 



states is calculated by recurrence relations of the 
efficient numerical algorithm motivated by Ref. 
quantum problems, open or closed. It could give results in a much shorter time compared 
with the traditional methods, such as the well-known 4-order Runge-Kutta algorithm, 
under the same requirement of numerical accuracy. 

After some derivations, the density matrix of the whole system p{t) can be determined 
by Eqs. [12] and [181 Tracing out the degrees of freedom of the environment, we finally obtain 
the dynamics of the subsystem qubits: 

p,(t)=Tr,(p(t)). (20) 
III. SIMULATION RESULTS AND DISCUSSIONS 



With Ps{t), we can discuss: (i) the concurrence 4l|, |42|, which is a very good measurement 



for the intra-entanglement of two two-level particles and defined as: 

C = max{Ai - A2 - A3 - A4, 0}, (21) 



where Aj, i = 1,2, 3, 4, are the square roots of the eigenvalues of the product matrix ps{a 



y , 



ay)p*{o'y (g) cr^) in decreasing order; (ii) the fidelity 43|, which is defined as 



Fd{t) = Tr[pideai(t)p.(t)]. (22) 

Pideai(^) represents the pure state evolution of the subsystem under Hg only, without 
interaction with the environment. The fidelity is a measure for decoherence and depends 
on Pideai- It achieves its maximum value 1 only if the time dependent density matrix ps{t) 



is equal to Pidcai(t)- 



The results and discussions about the two quantities are divided into two subsections, 
IIII Al and IIII B[ Both of these two cases suggest a type of controlling method to the 
entangled qubits, especially the latter one. 

It is interesting to note that no matter which one of the four Bell states is chosen as the 
initial state for the subsystem, there is no distinction between their dynamics of concurrence. 
So we need not to concretely point out the initial states in the following discussions. But 
there are differences between the fidelity dynamics of jei^a) and |e2,4). It is shown that the 
spins in the initial state of the 1st group, |ei) and 63), are parallel; while in that of the 
2nd group, |e2) and 64), are antiparallel. That is why the two groups have different fidelity 
evolution behaviors. In fact, any other basis obtained from the Bell states by replacing the 



coefficients ±1/^2 with e*7V2 {9 being a real number) has the same dynamics [4J]. 



A. Thermodynamic limit {N 00) 

The free evolution (It means the subsystem is decoupled from the bath or = 0.) 
of both concurrence and fidelity of the subsystem is time-independent, because (i) Hg 
in Eq. [T] {Hs = fJ'oi'^A + ^b)) "^^^ ^'^^ change the entanglement degree between the 
two qubits; (ii) when qq = 0, the effect of the bath is excluded from the evolution of 
the subsystem. Thus if go = 0, we always have C{t) = Fd{t) = 1. In Fig. 1(a), we 
show different effects on the dynamics by four anisotropic parameters: 7 = 0,0.2,0.6,1. 
From the four curves, we obtain two findings, (i) When 7 is not too large, the entan- 
glement (C(0) = 1) can always be recovered to a high degree after some oscillations. 
For instances, ^(7 = 0,got = 8.960) = 0.974631, C(7 = 0.2, = 8.912) = 0.961513, 
C(7 = 0.6, got = 12.672) = 0.844158. Yet the concurrence will never reach 1 in a long time 
scale, (ii) The curve of 7 = 1 shows a totally different behavior from the other three cases. 
It keeps decreasing with some little fluctuations. Then we turn to the other sub-flgures. 



Fig. 1(c) is almost the same as Fig. 1(b) There are some obvious disagreements between 



Fig. 1(b) and Fig. 1(a) For example, ^(7 = 0,got = 8.960) = 0.974631, but at the 



same time, Fdl{'-f = 0,got = 8.960) = 0.101803. In other word, when the concurrence 
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FIG. 1: Time evolution of (a) Concurrence from all 4 Bell states, (b) Fidelity from jei^a), (c) 
Fidelity from |e2,4) at different anisotropic parameter: 7 = (solid curve), 7 = 0.2 (dashed curve), 
7 = 0.6 (dot dashed curve), 7 = 1.0 (dotted curve). Other parameters are N ^ 00, fiQ = 2go, 
9 = 9o,T = go- 

of the subsystem has been mostly retrieved, the state of that is not simultaneously back 
to the its initial state. Only the combination of the concurrence and the fidelity can give 
a complete description of the real revival of the state. It can be testified in the case of 
7 = (solid curves in the two sub- figures) . When got = 15.624, the concurrence evolves to 
C = 0.875742 and the fidelity goes back to Fdl = 0.924277. So at that moment, Ps{got) is 
mainly composed by |ei,3)(ei,3|. 
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{c)Fd2{t) 

FIG. 2: Time evolution for (a) Concurrence from all 4 Bell states, (b) Fidelity from jei^s), (c) 
Fidelity from |e2,4) at different values of temperature: T = 0.2go (solid curve), T = go (dot dashed 
curve), T = 5go (dashed curve). Other parameters are — > oo, = 2^07 5 = 50; 7 = 0.6. 



In Fig. 2(a), Fig. 2(b) and Fig. 2(c), we plot the dynamics of the concurrence and 



fidelity at different temperatures. When the temperature is not too high, such as T = 0.2gQ 
and T = Igo, both concurrence and fidelity represent a periodical oscillation. At some 
moments, they can restore to a high degree. The restoring degree of the both quantities, 
however, decreases as the temperature increases. Similar to Fig. [T], the revivals of the 
concurrence and fidelity do not take place simultaneously. When the bath is at a high 
temperature, such as T = 5go, the concurrence quickly declines to zero (to see the dashed 



line in Fig. 2(a)) and does not go back to C > immediately. It means that when the 
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FIG. 3: Time evolution for (a) Concurrence from all 4 Bell states, (b) Fidelity from jei^s), (c) 
Fidelity from |e2,4) at different values of coupling strength between subsystem and bath: g = 2gQ 
(solid curve), g = 4(/o (dashed curve), g = 8go (dot dashed curve). Other parameters are N ^ oo, 
fJ-o = 2^0, 7 = 0-6, T = go. 



local spin bath is adjusted to a high temperature, it makes a sudden disappearance to the 
entanglement of a non-localized state and it will lose the control ability to the subsystem. 



45 



46j. In 



This is the effect that has been called "entanglement sudden death" (ESD) 
Ref. {3], after the concurrence goes abruptly to zero, it arises more or less from nowhere, 
since there is no local effect under the action of weak noises. Our model is still an ex- 



ample of ESD, however, the concurrence arises after some time due to the local thermal bath. 
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To find out the role of the subsystem-bath couphng strength g, we keep the bath at a 



moderate temperature T = g^. In Fig. 3(a), all the three curves show periodical behaviors 



and the oscillation amplitudes are strikingly damped by increasing g from 2gQ to 8gQ. When 
g is up to 8gQ, the fiuctuation magnitude of concurrence near C = 1 is too small to be 
noticed. It is like the case of go = 0, in which the bath is decoupled from the subsystem 
and C{t) = 1. It is consistent with the claims in Refs. 18,[l9|, [2^ that enough strong intra- 



coupling strength among bath spins can make the evolution of the subsystem be completely 
determined by the Hamiltonian of itself Hg. But the subsystem state does not receive the 



same protection as the subsystem entanglement degree, especially when g = Sgo. Fig. 3(b) 



manifests that the revival period of fidelity is much longer than that of concurrence. In 
fact, because Qubit B is not under the infiuence of the bath, the bath can not make a 
decoherence-suppression effect on the subsystem. 



B. finite bath spins (iV = 40) 



In the previous two-center-spin-spin-bath works [17|, |29|, it is supposed that the number 
of bath spins is infinite, which helps to reduce the Hamiltonian Eqs. [10] and [11] to a simple 
form. Yet as the controlling device in a real quantum information equipment, the spin 
bath, in principle, should be made of finite number of spins-1/2. Then in this subsection, 
we use the 1st order expansion of the Hamiltonian (Eqs. [TO]and [TTl) to introduce a finite N 
in the present problem. The error about this approximation is about 0(1/ N"^) as Eq. [11] 
indicates. Without loss of generality, we set = 40 . 



Comparing the result of Fig. [T] with that of Fig. |U we can find some agreements 
and some disagreements. The most identical characteristic between them is that the 
concurrence dynamics is independent of the choice of state as long as it is one of the 
four Bell states. Yet when 7 = 1, the concurrence (dotted curve) does not decrease 
monotonously during the given time. For fidelity, it is shown that with a bigger anisotropic 



parameter 7, the curves in Fig. 4(b) oscillate with a shorter period, which is opposite to 



the tendency in Fig. 1(b), While Fig. 4(c) is almost the same as Fig. 1(c), with a little 



longer oscillation period. It seems that the fidelity evolution of the 2nd group is not very 
sensitive to the bath-spin number N. The differences between the infinite and finite 
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FIG. 4: Time evolution for ((a) Concurrence from all 4 Bell states, (b) Fidelity from jei^s), (c) 
Fidelity from |e2,4) at different values of anisotropic parameter: 7 = (solid curve), 7 = 0.2 
(dashed curve), 7 = 0.6 (dot dashed curve), 7 = 1.0 (dotted curve). Other parameters are = 40, 
fJ-o = 2go, g = go,T = go. 

N cases might arise from their different energy-level numbers and corresponding weights 
in our numerical scheme (in subsection IIIIB|) . Under the same requirement of numeri- 
cal accuracy, for the infinite N, we need to consider 14, 15, 18 and 20 energy levels when 
7 is 0, 0.2, 0.6 and 1 respectively; for N = 40, we calculate 9, 10, 17 and 18 levels respectively. 



In the comparison of Fig. [5] with Fig. [2], we can also find the effect of a finite A^. At 
low temperature (T = 0.2go), the entanglement degree of the subsystem qubits oscillates 
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FIG. 5: Time evolution for (a) Concurrence from all 4 Bell states, (b) Fidelity from jei^s), (c) 
Fidelity from |e2,4) at different values of temperature: T = 0.2go (solid curve), T = go (dot dashed 
curve), T = 5go (dashed curve). Other parameters are = 40, fiQ = 2gQ, g = go, ^ = 0.6. 



with a nearly perfect period between the value of 0.8 and 1.0 (the solid curve in Fig. 



5(a)). However, the subsystem of Group 1 (Group 2) goes back to its own initial state 
only once in almost five (ten) revival periods of concurrence, which is illustrated by the 



corresponding curve in Fig. 5(b) (Fig. 5(c)). It is obvious that the increase of the 
temperature will also destroy this perfect oscillation. When T = Igo, the second peak 
value, C{got = 16.19) = 0.918798 is lower than the first one C{gQt = 0) = 1 (to see the 



dot dashed curve in Fig. 5(a)) and the peak Fdl{got = 11.55) = 0.877301 is higher than 



Fdl^got = 23.13) = 0.779953 (to see the dot dashed curve in 5(b)). When the temperature 
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(a)C(t) 



(b)Fdl(t) 




{c)Fd2{t) 

FIG. 6: Time evolution for (a) Concurrence from all 4 Bell states, (b) Fidelity from jei^s), (c) 
Fidelity from |e2,4) at different values of coupling strength between subsystem and bath: g = 2gQ 
(solid curve), g = Ago (dashed curve), g = 8go (dot dashed curve). Other parameters are = 40, 
^0 = 2^0, 7 = 0-6, T = go. 

is up to bgo, the entanglement vanishes to zero in a fairly short stretch of time. The entangle- 
ment "death" time is longer than that in the case of — oo (Fig. [2]). So whether N oo 
or N is finite, the revival of the concurrence after ESD results from the effect of thermal bath. 



For the subsystem-bath coupling g, the dynamics of concurrence in the case of = 40 



(to see Fig. 6(a)) is almost the same as that in the N ^ oo case (to compare it with Fig. 



3(a)). The evolution of the fidelity of the 1st group shows significant changes when N is 
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changed from infinity (Fig. 3(b) ) to 40 (Fig. 6(b) ) while the fidehty of the 2nd group does 
not show very obvious changes. And for the 1st group, all the three cases behave periodical 
oscillations. They manifest that the subsystem in the condition of = 40 can be restored 
to the initial state with more chances or possibilities than that in the condition of — oo. 



IV. CONCLUSION 



We studied the time evolution of two separated qubit spins with a thermal equilibrium 
bath composed of infinite or finite spins in a quantum anisotropic Heisenberg XY model. 
The bath can be treated effectively as a single pseudo-spin of N/2 according to the symmetry 
of the Hamiltonian. By the Holstein-Primakoff transformation and the first order of 1/A^ 
expansion, it is further considered as a single-mode boson field. The pair of qubits served 
as an quantum information device is initially prepared in a Bell state. It is interesting 
that the concurrence and the fidelity dynamics of the subsystem can be controlled by some 
characteristic parameters of the spin bath. Through the adjustment, we show that (i) the 
concurrence dynamics of the subsystem is independent of the initial state, whether A^ is 
infinite or finite, however, the fidelity dynamics is divided into two groups; (ii) smaller 
anisotropic parameter 7 can help the subsystem to evolve into a highly-entangled state, but 
this restoration should be measured by the combination of concurrence and fidelity; (iii) 
the bath at higher temperature makes a sudden death to the entanglement (ESD) of the 
subsystem and strongly destroys the fidelity of that; (iv) the spin-bath can help to keep 
the high entanglement degree between the two subsystem spins in the condition of large 
intra-coupling g. 
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